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Abstract 

The aim of this paper is two-fold. First we analyze the sequence of intensity measures 
of a spatial branching point process arising in a multiple target tracking context. We 
study its stability properties, characterize its long time behavior and provide a series of 
weak Lipschitz type functional contraction inequalities. Second we design and analyze 
an original particle scheme to approximate numerically these intensity measures. Under 
appropriate regularity conditions, we obtain uniform and non asymptotic estimates and 
a functional central limit theorem. To the best of our knowledge, these are the first 
sharp theoretical results available for this class of spatial branching point processes. 

Keywords : Spatial branching processes, multi-target tracking problems, mean field 
and interacting particle systems, w.r.t. time, functional central limit theorems. 

1 Introduction 

Multi-target tracking problems deal with tracking several targets simultaneously given noisy 
sensor measurements. Over recent years, point processes approaches to address these prob- 
lems have become very popular. The use of point processes in a multiple-target tracking 
context was first proposed in S. Mori et al. |15] as early as in 1986. Using a random sets 
formalism, a formalism essentially equivalent to the point process formalism fl6], R. Mahler 
and his co-authors proposed in two books [11^ [12] a systematic treatment of multi-sensor 
multi-target filtering problems. However, as mentioned in |16] . "... although the random 
sets formalism (or the point process formalism) for multitarget tracking has provided a uni- 
fied view on the subject of multiple target tracking, it has failed to produce any significant 
practical tracking algorithms...". 

This situation has recently changed following the introduction of the PHD (probability 
hypothesis density) filter by R. Mahler [T^ . The PHD filter is a powerful multi-target 
tracking algorithm which is essentially a Poisson type approximation to the optimal multi- 
target filter [TSldlllT?]. It has found numerous applications since its introduction. The 
PHD filter cannot be computed analytically but it can be approximated by a mixture of 
Gaussians for linear Gaussian target models |18) and by non-standard particle methods for 
nonlinear non-Gaussian target models [9l llOj. 
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Despite their increasing popularity, the theoretical performance of these multi-target 
particle methods remain poorly understood. Indeed their mathematical structure is signifi- 
cantly different from standard particle filters so the detailed theoretical results for particle 
filters provided in [3] are not applicable. Some convergence results have already been es- 
tablished in [9l [To] but remain quite limited. Reference [9] presents a basic convergence 
result for the PHD filter but does not establish any rate of convergence. In ^Oj the authors 
provide some quantitative bounds and a central limit theorem. However these quantitative 
bounds are not sharp and no stability result is provided. 

The aim of this work is to initiate a thorough theoretical study of these non-standard 
particle methods by first characterizing the stability properties of the "signal" process and 
establishing uniform w.r.t. the time index convergence results for its particle approximation. 
This "signal" process is a spatial branching point process whose intensity measure always 
satisfies a closed recursive equation in the space of bounded positive measures. We will not 
consider any observation process in this article. The analysis of the particle approximations 
of PHD filters is presented in [7j. It builds heavily upon the present work but is even more 
complex as it additionally involves at each time step a nonlinear update of the intensity 
measure. 

The rest of this paper is organized as follows: 

In section [21 we present a spatial branching point process which is general enough to 
model a wide variety of multiple target problems. We establish the linear evolution equation 
associated to the intensity measures of this process and introduce an original particle scheme 
to approximate them numerically. Section [3] summarizes the main results of this paper. 
In Section [H we provide a detailed analysis of the stability properties and the long time 
behavior of this sequence of intensity measures, including the asymptotic behavior of the 
total mass process, i.e. the integral of the intensity measure over the state space, and the 
convergence to equilibrium of the corresponding sequence of normalized intensity measures. 
For time-homogeneous models, we exhibit three different types of asymptotic behavior. The 
analysis of these stability properties is essential in order to guarantee the robustness of the 
model and to obtain reliable numerical approximation schemes. Section [5] is devoted to 
the theoretical study of the non-standard particle scheme introduced to approximate the 
intensity measures. Our main result in this section is a non-asymptotic convergence for 
this scheme. Under some appropriate stability conditions, we additionally obtain uniform 
estimates w.r.t. the time parameter. 

2 Spatial branching point process and its particle approxi- 
mation 

2.1 Spatial branching point process for multi-target tracking 

Assume that at a given time n there are A^^^ target states (A'^)i<i<Af„ taking values in some 
measurable state space En enlarged with an auxiliary cemetery point c. The state space En 
depends on the problem at hand. It may vary with the time parameter and can include all 
the characteristics of a target such as its type, its kinetic parameters as well as its complete 
path from the origin. As usual, we extend the measures 7^ and the bounded measurable 
functions fn on En by setting 7ra(c) = and /n(c) = 0. 

Each target has a survival probability e„(A'^) G [0,1]. When a target dies, it goes to 
the cemetery point c. We also use the convention en(c) = so that a dead target can only 
stay in the cemetery. Survival targets give birth to a random strictly positive number of 
individuals /i^(X^) where (/i^(X^)) is a collection of independent random variables 

such that E (/i^(x„)) = Hn{xn) for any x„ E En where Hn is a given collection of bounded 
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functions H^- We have Hn (xn) > 1 for any x„ £ En as h^ixn) > 1- This branching 
transition is called spawning in the multi-target tracking literature. We define G„ = e„if„. 

After this branching transition, the system consists of a random number Nn of individuals 
. Each of them evolves randomly = Xn ^n+i according to a Markov 
transiFion Mn+i{xn,dxn+i) from En into En+i- We use the convention Mn+i{c^c) = 1, so 
that any dead target remains in the cemetery state. 

At the same time, an independent collection of new targets is added to the current 
configuration. This additional point process is modeled by a spatial Poisson process with 
a prescribed intensity measure on En+i- It is used to model new targets entering the 
state space. 

At the end of this transition, we obtain A^n+i = Mi + ^'n+i targets (X^_|_i)i<j<Ar„_,_i, 
where A^^+i is a Poisson random variable with parameter given by the total mass //n+i(l) of 

the positive measure /in+i, and (-'^^i^*)i<i<Ar'^^ are independent and identically distributed 
random variables with common distribution 7i„-|_i = ^n+i//Jn+i(l) where := fin+i{dx). 

Example. To illustrate the model, we present here a simple yet standard example JT^ 
of a target evolving in a two-dimensional surveillance region S CM?. In this case, we set 
En = E = S xM?. All the targets are assumed to be of the same type. The state of a target 
-^n = \Pn^Pn,Vn,Vn]^ consists of its position {Pn,Pn) S S and velocity {Vn,Vn) G M? and is 
assumed to evolve according to a linear Gaussian model 

Xn = AXn-1 + Vn (2.1) 

where A is a known transition matrix and Vn ~ A/'(0, S) is a sequence of i.i.d zero-mean 
normal random variables of covariance T,; i.e. Mn{xn-i,dxn) = M{xn-i,Xn)dxn with 

M{Xn-l,Xn) = |27rSp^/^exp ^-^ {Xn - AXn-l)^^~^ (x„ - AXn-l)^ ■ 

In the example, we assume that fj.n {x) = fi{x), e„(a;) = s > and hn{xn) = /i G {1, 2} with 
P (/i = 1) = 1 — P (/i = 2) = a. Hence for this model, each target Xn-i survives at time 
n — 1 with a probability s. Each survival target has one offspring with probability a which 
evolves according to 1^2. 1\) or two offspring with probability \ — a which, conditional upon 
Xn-i, independently evolve according to 112.1]) . Additionally, a random number of targets 
distributed according to a Poisson distribution of parameter /u(l) appear. These targets are 
independent and distributed in E asji = fi/ fi{l). 



2.2 Sequence of intensity distributions 

At every time n, the intensity measure of the point process Xn := Yli^=i ^x^ associated to 
the targets is given for any bounded measurable function / on £"„ U {c} by the following 
formula: 

7„(/) :=E(A'„(/)) with Xn{f):=l f{x)Xn{dx) 

To simplify the presentation, we suppose that the initial configuration of the targets is a 
spatial Poisson process with intensity measure fiQ on the state space Eq. 

Given the construction defined in section 12.11 it follows straightforwardly that the inten- 
sity measures 7„ on En satisfy the following recursive equation. 

Lemma 2.1 For any n >0, we have 

jn+i{dx') = -fn{dx) Qn+i{x,dx') -\- ^n+i{dx') (2.2) 
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with the initial condition 70 = /Uq where ^n+i is the intensity measure of the spatial point 
process associated to the birth of new targets at time n + 1 while the integral operator Qn+i 
from En into En+i is defined by 



Qn+l{Xnj dXn+l) '■— Gn{Xn) M„+i ) . 

Proof: 

For any bounded measurable function / on En+i U {c}, we have 



7„+i (/) = E [ / + E ^ / 

where, thanks to the Poisson assumption, we have 



E 



V 



and 



i=l 



•F n 



= e|e| J^M„+i(/)(x 



E ^ e„ (X;) i/„(X;)M„+i (/) {XI, 



(2.3) 



= 7n (e„i?„M„+i (/)) 

where J>i denotes the cr-field generated by {X^)-^^.^j^ and Qn the ci-field generated by 



These intensity measures typically do not admit any closed-form expression. A natural 
way to approximate them numerically is to use a particle interpretation of the associated 
sequence of probability distributions given by 

r]n{dx) := jn{dxn)/7n{'i^) with 7„(1) := / jn{dx) 

J En 

To avoid unnecessary technical details, we further assume that the potential functions G„ 
are chosen so that for any x G E^ 



< Qn - < Gn{x) < gn,+ < OO 



(2.4) 



for any time parameter n > 0. Note that this assumption is satisfied in most realistic multi- 
target scenarios such as the example discussed at the end of section [2TTJ Indeed the condition 
On - < Gn{x) essentially states that there exists e„^_ > such that e„ (x) > e^,- for any 
x ^ En as Hn {x) > 1. The condition Gn{x) < gn,+ states that there exists -ffn,+ < 00 such 
that Hn (x) < Hn,+ for any x G En as e„ (x) < 1. In the unlikely scenario where (|2.4p is 



4 



not satisfied then the forthcoming analysis can be extended to more general models using 
the techniques developed in section 4.4 in see also [3]. We denote by V{En) the set of 
probability measures on the state space E^- 

To describe these particle approximations, it is important to observe that the pair process 
(7„(l),7y„) e (M+ X V{En)) satisfies an evolution equation of the following form 

(7n(l),?/n) = r„(7n-l(l),r/n-l) (2.5) 

We let and be the first and the second component mappings from (M-|_ x V{En)) into 
M+, and from (M_|_ x V{En)) into V{En)- The mean field particle approximation associated 
with the equation (j2.5p relies on the fact that it is possible to rewrite the mapping T'^^^ in 
the following form 

r^+l(7n(l),r/n) = ?/ni^n+l,(7„(l),r,„) (2-6) 

where i^n+i,(m,r?) is a Markov kernel indexed by the time parameter n, a mass parameter 
m G M_|_ and a probability measure rj on the space En- In the literature on mean field particle 
systems, Kn,(m,ri) is called a McKean transition. The choice of such Markov transitions 
Kn,(m.-q) is not Unique and will be discussed in section [5Tl 
Before concluding this section, we note that 

7„+i(drc') = {'^nQn+i){dx') := / -fn{dx) Qn+i{x,dx') (2.7) 



when = 0. In this particular situation, the solution of the equation (12. 2p is given by the 
following Feynman-Kac path integral formulae 

7n(/) = 7o(l) E I /(^n) n G^p)] (2-8) 

y 0<p<n j 

where Xn stands for a Markov chain taking values in the state spaces E^ with initial dis- 
tribution ryo = 70/70(1) and Markov transitions M„ (see for instance section 1.4. 4. in [4j). 
These measure- valued equations have been studied at length in P]. 

2.3 Mean field particle interpretation 



The transport formula presented in (12. 6j) provides a natural interpretation of the probability 
distributions rjn as the laws of a process X„ whose elementary transitions Xn -'^n+i de- 
pends on the distribution r/„ = Law(X„) as well as on the current mass 7n(l)- In contrast to 
the more traditional McKean type nonlinear Markov chains presented in [3| , the dependency 
on the mass process induces a dependency of the whole sequence of measures r^p, from the 
origin p = up to the current time p = n. 

From now on, we will always assume that the mappings 



1<^<NJ ^ ^ ^n) ^ ^^^1^5^,) {x^An+l] 



are measurable w.r.t. the product sigma fields on (M_(_ x E^), for any > 0, > 1, and 
1 < i < A^, and any measurable subset An+i C En+i- In this situation, the mean field 
particle interpretation of (j2.6p is an E^ -valued sequence ei^) = f^^^'^)) defined as 



l<i<Af 



7^+1(1) = 7^(l)^^(G„) + Mn+l(l) 



ln+1 

(2.9) 

[(.i+l^dx = nili ^n+i,(7^(i),,,^)(d^'*\rf2;*) 
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with the pair of occupation measures {jn^Vn) defined below 

1 ^ 

nYI '^ei'^-') '^n {dx) := 7^ (1) 7]^ (dx) 

1=1 

In the above displayed formula, J-n^^ stands for the u-field generated by the random sequence 
{£,p^^)o<p<n, and dx = dx^ x . . . x dx^ stands for an infinitesimal neighborhood of a point 
X = (x^, . . . ,x^) G . The initial system ^q^'' consists of N independent and identically 
distributed random variables with common law ijq. As usual, to simplify the presentation, 
we will suppress the parameter when there is no possible confusion, so that we write 
and instead of and d"^'^). 

In the above discussion, we have implicitly assumed that the quantities /^n(l) are known 
and that it is easy to sample from the probability distribution 7I„(dx) := /i„((ix)/;U„(l). 
In practice, we often need to resort to an additional approximation scheme to approximate 
and Jin- This situation is discussed in section [6l This additional level of approximation 
has essentially a minimal impact on the properties of the particle approximation scheme 
which can be analyzed using the same tools. 



2.4 Notation 

For the convenience of the reader, we end this introduction with some notation used in the 
present article. We denote by Ai{E) the set of measures on some measurable state space 
{E,£) and we recall that ViE) is the set of probability measures. We also denote 13{E) the 
Banach space of all bounded and measurable functions / equipped with the uniform norm 
11/11 and Osci(i?) the convex set of £^-measurable functions / with oscillations osc(/) < 1 
where osc(/) = sup \f (x) - f {y)\. 

{x,y)&E'^ 

We let //(/) = / fi{dx) f{x) be the Lebesgue integral of a function / G S{E) with 
respect to a measure /x G Ai(E). We recall that a bounded integral kernel M{x, dy) from a 
measurable space {E, £) into an auxiliary measurable space {E\ £') is an operator / i— )• M{f) 
from B{E') into 13(E) such that the functions x i— )■ M{f){x) := J^, M{x, dy)f{y) are im- 
measurable and bounded for any / G B{E'). The kernel M also generates a dual operator 
/i ^ fj,M from M{E) into MiE') defined by {fJ.M){f) := fi{M{f)). A Markov kernel is 
a positive and bounded integral operator M with M(l) (x) = 1 for any x £ E. Given a 
pair of bounded integral operators (Mi,M2), we let (M1M2) be the composition operator 
defined by (MiM2)(/) = Mi(M2(/)). For time-homogenous state spaces, we denote by 
j^k ^ M'^-^M = MM^-^ the fc-th composition of a given bounded integral operator M, 
with A; > 0, with the convention = Id the identity operator. We also use the notation 

M ([/i - M(/i)] [h - M(/2)]) (x) := M ([/i - M(/i)(x)] [/2 - M(/2)(x)]) (x) 

for some bounded functions /i, /2. 

We also denote the total variation norm on M.{E) by ||/i||tv = sup^gQ^^^^^;) When 
the bounded integral operator M has a constant mass, that is M(l) (x) = M(l) (y) for any 
(x, y) G -B^, the operator n 1— )• maps M.{E) into M.{E'). In this situation, we let P{M) 
be the Dobrushin coefficient of a bounded integral operator M defined by the following 
formula 

/3(M) := sup {osc(M(/)) ; / G Osci(i?)} 

Given a positive function G on E, we let ■ V ^ ^(-^) ^ ^civ) ^ 'PiE) be the 

Boltzmann-Gibbs transformation defined by 

^'G(r?)(dx) := — ^ G{x) r]{dx) 
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We recall that ^civ) can be expressed in terms of a Markov transport equation 

TjSrj = ^g{v) (2.10) 

for some selection type transition Srj{x,dy). For instance, for any e > s.t. G{x) > e for 
any x, we notice that 

-Jfi -^iy) 



so we can take 



S^{x,dy) := ^ S,{dy) + (^1 - ) ^ ^G-e)iv){dy) (2.11) 



For e = 0, we have S.^{x,dy) = '$Q{r]){dy). We can also choose 

S^{x, dy) := eG{x) 6,{dy) + (1 - eG{x)) ^G{v){dy) (2.12) 

for any e > that may depend on the current measure t], and s.t. eG{x) < 1. For instance, 
we can choose 1/e to be the r^-essential supremum of G. 

3 Statement of the main results 

At the end of section 12.21 we have seen that the evolution equation (j2.2p coincides with that 
of a Feynman-Kac model (12. 8p for /x„ = 0. In this specific situation, the distributions 7„ 
are simply given by the recursive equation 

7„ = 7„_i(5„ =^ VO < p < n 7„ = 7p(5p,„ with Qp^n = Qp+i ■ ■ ■ Qn-iQn (3.1) 

For p = n,we use the convention Qn,n = Id. In addition, the nonlinear semigroup associated 
to this sequence of distributions is given by 

Vnif) = '^p,n{Vp){f) ■= VpQpAf)/VpQp,n{'^) = Vp {Qp,n{'^)Pv,nU)) MpQp^niX) (3.2) 

with the Markov kernel Pp^ni^p, dxn) = Qp,n{xp,dxn) / Qp,n{xp, En). The analysis of the 
mean field particle interpretations of such models has been studied in [4j. Various prop- 
erties including contraction inequalities, fluctuations, large deviations and concentration 
properties have been developed for this class of models. In this context, the fluctuations 
properties as well as L^-mean error estimates, including uniform estimates w.r.t. the time 
parameter are often expressed in terms of two central parameters: 

qp^n = sup ^^'"/lll^l and /3(Pp,„) = sup ||Pp,„(x, .) - ^'p,n(y, Olltv (3.3) 

x,y '^p,n{i')(y ) x,y€Ep 

with the pair of Feynman-Kac semigroups {Pp,n,Qp,n) introduced in (j3.ip and (j3.2p . 
We also consider the pair of parameters {g-{n), g+{n)) defined below 

g-{n) = inf inf Gp < sup s\xpGp = g+{n) 

0<p<n Ep 0<p<n Bp 

We also write g_/^{n) to refer to both parameters. The first main objective of this article is 
to extend some of these properties to models where /i„ is non necessarily null. We illustrate 
our estimates in three typical scenarios 

1) G = 5-/+ = 1 2) 5+ < 1 and 3) g- > 1 (3.4) 
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arising in time homogeneous models 

{En,Gn, Mn, Hn, g~{n), g+{n)) = {E,G,M,fi,g^,g+) (3.5) 

These three scenarios correspond to the case where, independently from the additional spon- 
taneous births, the existing targets die or survive and spawn in such a way that either their 
number remains constant (G = g-/+ = 1), decreases (5+ < 1) or increases {g- > 1). 

Our first main resuh concerns three different types of long time behavior for these three 
types of models. This result can basically be stated as follows. 

Theorem 3.1 For time homogeneous models i3. 5\) . the limiting behavior 0/ (7„(1), ry^) in 
the three scenarios is as follows: 

1. When G{x) = 1 for any x £ E, we have 

7„(1) = 7o(l) + ^(1) n and - r/oo||tv = O 

when M is chosen so that 

sup ||M"(x, .) — ??oo||tv < 00 for some invariant measure r/00 = "HooM. (3.6) 

2. When g+ < 1, there exists a constant c < 00 such that 

V/ G B{E), |7„(/) - 7oo(/)| V \vnif) - r/oo(/)| < c 5+ 11/11 
with the limiting measures 

7oo(/) := j;^Q"(/) and i^^f) := 7oo(/)/7oo(l) (3.7) 

n>0 

3. When g- > 1 and there exist k > 1 and e > such that M^{x, .) > e M^{y, .) for any 
x,y £ E then the mapping $ = ^n-i,n introduced in f3.^) has a unique fixed point 
Voo = ^(^00) and 

lim ilog7„(l) = logr?oo(G') and ||r/„ - 7?oo||tv < c e~^" 

n— >oo n 

for some finite constant c < 00 and some A > 0. 

A more precise statement and a detailed proof of the above theorem can be found in 
section I4.2[ 

Our second main result concerns the convergence of the mean field particle approxi- 
mations presented in (j2.9p . We provide rather sharp non asymptotic estimates including 
uniform convergence results w.r.t. the time parameter. Our results can be basically stated 
as follows. 

Theorem 3.2 For any n > 0, and any N > 1, we have 7n(l) and 7^ (1) G I„ with the 
compact interval In defined below 

n 

I„ := [m_(n),m+(n)] where m^/^{n) := /Xp(l)g_/+(re)("~P) (3.8) 

p=0 
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In addition, for any r > 1, f £ Osci(-E„), and any N > 1, we have 
VNE(\[rj^-r^r.] if) 



1 

I r \ r 



< Or bn with bn < 6p,^ 
p=0 



(3.9) 



where Or < oo stands for a constant whose value only depends on the parameter r and bp^n 
is the collection of constants given by 



bp^n '■ — 2 (1 A TTlp^n) Qp,n 



with the pair of parameters 



p<q<n ^P<r<n ^r,n 



q,nj 



(3.10) 



m, 



= m+(p)||(5p,„(l)||/ ^ Cg^n and Cp^n l^pQpA^) 

p<q<n 



Furthermore, the particle measures 7^ are unbiased, and for the three scenarios \3.4^ 
with time homogenous models s.t. M^{x, •) > e M^[y, .), for any x,y £ E and some pair of 
parameters k > 1 and e > 0, the constant bn in ^3. 9\) can be chosen so that su.p^>Q bfi < 00/ 
in addition, we have the non asymptotic variance estimates for some d < 00, any n > 1 and 
for any N > 1 



E 



7^ 

7n(l) 



< d 



n+1 
N -1 



1 + 



d 



N -I 



(3.11) 



The non asymptotic estimates stated in the above theorem extend the one presented 
in [3l H] for Feynman-Kac type models p.Sp where fin = 0. For such models, the L^-mean 
error estimates (|3.9p are satisfied with the collection of parameters bp^n ■= "^Qp^n PiPp,n), 
with p < n. The extra terms in (I3.10p are intimately related to /i„ whose effects in the 
semigroup stability depend on the nature of Gn ■ We refer to theorem 13. H section 14.21 and 
section 14.31 for a discussion on three different behaviors in the three cases presented in ()3.4p . 

A direct consequence of this theorem is that it implies the almost sure convergence 
results: 



lim < (/) = 77„(/) and lim 7/,^ (/) = 7„(/) 

for any bounded function / G B{En)- 

Our last main result is a functional central limit theorem. We let W:^ be the centered 
random fields defined by the following formula 



N 

^n-l^n,(7,^'(l),<„i) 

We also consider the pair of random fields 



+ 



1 



N 



N 



(3.12) 



N 



iV[7?^-r?„] and V^'"' := 



N 



For n = 0, we use the convention W, 



N 



Theorem 3.3 The sequence of random fields {Wn)n>o converges in law, as N tends to 
infinity, to the sequence of n independent, Gaussian and centered random fields iyVn)n>o 
with a covariance function given for any f,gG B{En) and n >0 by 



E{Wn{f)Wnig)) 

= ??n-li^n,(7u-i(l),)?n-i) ([/ - -f^n,(7n-i(l),r;„-i)(/)][9 " ^n,(7„_i(l),7;„_i) (s)]) ) • 



(3.13) 
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In addition, the pair of random fields Vn' and Vn converge in law as N ^ oo to a pair 
of centered Gaussian fields Vn and Vn defined by 

n , 
p=0 V7nUJ 

The details of the proof of theorem 13.21 and theorem l3.3l can be found in section [5^21 The 
proof of the non-asymptotic variance estimate (j3.1ip is given in section 15.2.11 dedicated to 
the convergence of the unnormalized particle measures 7^. The L^-mean error estimates 
()3.9p and the fluctuation theorem l3.3l are proved in section [5. 2. 2 i Under additional regularity 
conditions, we conjecture that it is possible to obtain uniform estimates for theorem 13.31 but 
have not estabhshed it here. 

The rest of the article is organized as follows. 

In sectional we analyze the semigroup properties of the total mass process 7n(l) and the 
sequence of probability distributions r/„. This section is mainly concerned with the proof of 
theorem 13. 1[ The long time behavior of the total mass process is discussed in section 14.11 
while the asymptotic behavior of the probability distributions is discussed in section 14.21 
In section 14.31 we develop a series of Lipschitz type functional inequalities for uniform 
estimates w.r.t. the time parameter for the particle approximation. In section O we present 
the McKean models associated to the sequence (7„(l),r7n) and their mean field particle 
interpretations. Section 15.21 is concerned with the convergence analysis of these particle 
approximations. In section [5.2.1l we discuss the convergence of the approximations of 7„(1), 
including their unbiasedness property and the non asymptotic variance estimates presented 
in (|3.1ip . The proof of the L^-mean error estimates (j3.9p is presented in section I5.2.2[ The 
proof of the functional central limit theorem 13.31 is a more or less direct consequence of the 
decomposition formulae presented in section [5^2] and is just sketched at the end of this very 
section. 



4 Semigroup analysis 

The purpose of this section is to analyze the semigroup properties of the intensity measure 
recursion ()2.2p . We establish a framework for the analysis of the long time behavior of these 
measures and their particle approximations (j2.9p . First, we briefly recall some estimate of the 
quantities (qp^n, PiPp,n)) in terms of the potential functions G„ and the Markov transitions 
Mn- Further details on this subject can be found in [3], and in references therein. 

We assume here that the following condition is satisfied for some A; > 1, some collection 
of numbers G (0, 1) 

{M)k Mp^p+k{xp, .) > ep Mp^p^kiVp, ■) with Mp^p+fe = Mp+iMp+2 • • • Mp+fc (4.1) 

for any time parameter p and any pair of states {xp,yp) ^ Ep. It is well known that the 
mixing type condition (M)^ is satisfied for any aperiodic and irreducible Markov chains on 
finite spaces, as well as for bi-Laplace exponential transitions associated with a bounded drift 
function and for Gaussian transitions with a mean drift function that is constant outside 
some compact domain. We introduce the following quantities 

6p,n ■■= sup Yl {Gq{xg)/Gq{yq)) and := (5p+i,p+fc (4.2) 

p<q<n 

where the supremum is taken over all admissible pair of paths with transitions Mq where an 
admissible path (a;p_i, Xp+i, x„,_i) is such that nn<g<n-^ {xq^i,dxq) > 0. Under the 
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above conditions, we have 0| p. 140] 

[n/fej-l 



l3{Pp,P+n) < n {'^~4+ik/i+ik) and qp,p+n < Sp^p+k/^p (4.3) 



For time-homogeneous Feynman-Kac models we set e := and 5k ■= ^o.fci for any k > 0. 
Using this notation, the above estimates reduce to [U p. 142] 

qp,p+n < 6k/e and /3(Pp,p+n) < (l - e^Sk-i) ^"^""^ (4.4) 
4.1 Description of the models 

The next proposition gives a Markov transport formulation of r„ introduced in ()2.5p . 
Proposition 4.1 For any n> 0, we have the recursive formula 

7„+l(l) = 7n(l) ^n(G'„) + /U„+i(l) 

= ^'g„(??„)M„+i^(^^(i)^^^) 



(4.5) 



with the collection of Markov transitions Mn+i^(^m,ri) indexed by the parameters m G M+ and 
the probability measures r/ G V{En) given below 

Mn+i,{m,r,){x,dy) := an{m,r]) Mn+i{x,dy) + {1 - an{m,r])) Jln+iidy) (4.6) 

with the collection of [0, l]-parameters an {m^if} defined below 

mr]{Gn) 

an (m, rj) = -- 

mri(Gn) + 

Proof: 

Observe that for any function / G B{En+i), we have that 

^ 7n{GnMn+l{f)) + fln+lif) ^ Injl) Vn{GnMn+l{f)) + fin+l{f) 
""^^ ln{Gn) + /^n.+l(l) 7n(l) Vn{Gn) + Mn+l(l) 

from which we find that 

r]n+i = an (7n(l)>^?n) ^>„+i(r?„) + (1 - a„ (7„(1),77„)) 
From these observations, we prove (|4.5p . This ends the proof of the proposition. ■ 

We let Tn+i be the mapping from M+ x V{En) into M+ x V{En+i) given by 

r„+i(m,r/) = (r^+i(m,r/),r2^i(m,r/)) (4.7) 
with the pair of transformations: 

= m r/(G„) + ^„+i(l) and r^+i(m,r/) = '^Gn{v)Mn+l,{m,7^) 
We also denote by (rp,n)o<p<„ the corresponding semigroup defined by 

\/0 < p < n Tp^n = ^p+l,vTp+l = ^n^n-l ■ ■ ■ Tp+l 

with the convention r„_„ = Id. 

The following lemma collects some important properties of the sequence of intensity 
measures 7„. 
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Lemma 4.2 For any < p < n, we have the semigroup decomposition 

In = lpQp,n + ^ /^gQg,n and 7„ = ^ lJLpQp,n (4.8) 
p<q<n 0<p<n 

In addition, we also have the following formula 

n 

7„(i) = /ip(i) n ^^(^^) (4-9) 

p=0 p<q<n 

Proof: 

The first pair of formulae are easily proved using a simple induction, and recalling that 
7o = /^o- To prove the last assertion, we use an induction on the parameter n > 0. The 
result is obvious for n = 0. We also have by 



7„+l(l) = 7„(5„+i(l) + = IniGn) + /in+l(l) 

This implies 

7n+l(l) = 7n(l) rin{Gn) + /^n+l(l) 

= 7„-l(l) ?/n-l(G„_i) r]n{Gn) + ^n(l) ??n(G'„) + Mn+l(l) 

n n+1 

=7o(i) n ^p(^p) +E '"p(i) n 

p=0 p=l p<g<n 

Recalling that jQ{dxo) = fj,o{dxo), we prove (j4.9p . This ends the proof of the lemma. ■ 

Using lemma 14. 2^ one proves that the semigroup Fp^^ satisfies the pair of formulae de- 
scribed below 

Proposition 4.3 For any < p < n, we have 

rp,„K??) =m r?Qp,„(l)+ ^ flqQq,ni'^) (4.10) 

p<q<n 

rp,„("^, r]) = Op^n {m, rf) $p,„(?7) + (1 - Qp.n ("i, 

(4.11) 

wit/i i/ie collection of parameters Cp^n ■= fJ'pQp,n{^) and the [0, l]-valued parameters Op^n ("^)^) 
defined below 



m 



Qp,n(-1) 



X]p<(j<ra '-ail 



(4.12) 



One central question in the theory of spatial branching point processes is the long time 
behavior of the total mass process 7n(l). Notice that 7n(l) = E(A'„(1)) is the expected size 
of the n-th generation. For time homogeneous models with null spontaneous branching /i„ = 
/i = 0, the exponential growth of these quantities are related to the logarithmic Lyapunov 
exponents of the semigroup Qp,n. The prototype of these models is the Galton- Watson 
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branching process. In this context three typical situations may occur: 1) 7^(1) remains 
constant and equals to the initial mean number of individuals. 2) 7n(l) goes exponentially 
fast to 0, 3) 7„(1) grows exponentially fast to infinity, 

The analysis of spatial branching point processes with /i„ = /i 7^ considered here is 
more involved. Loosely speaking, in the first situation discussed above the total mass process 
is generally strictly increasing; while in the second situation the additional mass injected in 
the system stabilizes the total mass process. Before giving further details, by lemma [42] we 
observe 7ri,(l) G In, foi' any n > 0, with the compact interval /„ defined in 13.81 

We end this section with a more precise analysis of the effect of // in the three scenarios 

(1321). 

In the further developments of this section, we illustrate the stability properties of the 
sequence of probability distributions ry„ in these three scenarios. 

1. When G{x) = 1 for any x £ E, the total mass process 7n(l) grows linearly w.r.t. the 
time parameter and we have 



7„(1) = m_(n) = m+(ri) = 70(1) + n 
Note that the estimates in ()4.12p take the following form 

7o(l)+;u(l)p 



(4.13) 



ap,n (7p(l),^p) < ap,„(7p(l)) := 1 A 



/x(l) (n-p) 







2. When < 1, the total mass process 7ra(l) is uniformly bounded w.r.t. the time 
parameter. More precisely, we have that 



V+W=5!!/+7o(l) + (l-5V) Y 



9-1+ 



This yields the rather crude estimates 

Ml 



7o(l) A 



1-5- 



< 7n(l) < 7o(l) V 



/^(l) 

l-5^ 



(4.14) 



We end this discussion with an estimate of the parameter Qp^„(m) given in ()4.12p . 
When the mixing condition (M)^ stated in ()4.ip is satisfied for some k and some fixed 
parameters Ep = e, using (j4.4p we prove that 



E 

p<r<n 



MQr,n(l) ^ 6^(1) 



'3p,r(Qr,n(l)) 



> 



E 



1 



p^n ^P''^^^) 



> 



e/^(l) 9 



-(n-p) 



1-3+ 



from which we conclude that for any n > p and any m £ Ip 
aljm) < 1 A 



(n-p) h (1 - 9+) 
m g _^_ 



eMi)(i-3f"^^) 



< 1 A 

< 1 A 



7o(l) V gt'^ 5,/{ef,{l)) 

i — .94- / 



*'(n— p)— )-oo 







(4.15) 



3. When g- > 1, the total mass process 7n(l) grows exponentially fast w.r.t. the time 
parameter and we can easily show that 

1 



5_ > 1 ^ 7„(1) > m_(n) = 7o(l) g^l + ^(1) 



9- 



1 



(4.16) 
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4.2 Asymptotic properties 

This section is concerned with the long time behavior of the semigroups Tp^n in the three 
scenarios discussed in (j4.13p , (j4.14p , and (j4.16p . Our results are summarized in theorem 13.11 
We consider time-homogeneous models (£'n, G^, M„, = {E,G,M,iJ,). 

1. When G{x) = 1 for any x £ E, we have seen in M.lSh that 7n(l) = 7o(l) + 
n. In this particular situation, the time-inhomogeneous Markov transitions 
,(7„_i(i),r)„_i) mtroduced m (|45]) are given by 

Mn{x,dy) =(l- —-i^ilL-^ M{x,dy) + Wv) 



7o(l) + nMl)y ' ' 7o(l) + nMl) 

This shows that r/„ = Law(X„) can be interpreted as the distribution of the states X„ 
of a time inhomogeneous Markov chain with transitions M„ and initial distribution 
If we choose in (j2.6p -ft^n+i,(7„(i),r?„) = ^n+i, the A^-particle model (j2.9p reduces 
to a series of independent copies of Xn- In this situation, the mapping Tq is given 
by 



7o(l) + n^(l) 7o(l)+n^(l) 

The above formula shows that for a large time horizon n, the normalized distribution 
flow rjn is almost equal to ^ X]o<p<n assume that the Markov kernel M is 
chosen so that (j3.6p is satisfied for some invariant measure r^oo = VooM. In this case, 
for any starting measure 70, we have 

II II ^ 70(1) nn{l) 1 ^ ^/l 

7o(l)+n^(l) 7o(l)+n^(l) n \n 

with T„ = sup^g£; ||M"(x, .)~^oo||tv For instance, suppose the mixing condition (M)^ 
presented in (14. ip is met for some k > 1 and e > 0. In this case, the above upper 
bound is satisfied with r„ = (1 — e)L"/'^J . 

2. Consider the case where 5+ < 1. In this situation, the pair of measures (13.7P are well 
defined. Furthermore, for any / € B{E) with ||/|| < 1, we have the estimates 

|7n(/) - 7oo(/)| < 7o(l) r?og"(l) + E -"^"(1) 

< [70(1) + Ai(l)/(1 - 9+)] -^n-.oo 
In addition, using the fact that 7n(l) > l^i^)^ we find that for any / G Osci(S) 

\Vn{f) - Voo{f)\ < -^TTT hn[f - Voo{f)] - 7oo[/ " Voo{f)]\ 
7n(l) 

< gl [7o(l)/;u(l) + 1/(1 - g+)] ^n^oo 

3. Consider the case where g- > 1. We further assume that the mixing condition (M)^ 
presented in (I4.ip is met for some k > 1 and some fixed parameters = e > 0. In 
this situation, it is well known that the mapping <I> = ^n-i,n introduced in (j3.2p has 
a unique fixed point r/00 = <^(r/oo), and for any initial distribution r]Q, we have 

W^oAvo) - VooWtv < a e-^ (4.17) 
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with 



as well as 



A = --log (l - e^/5o,fc-i) and a = 1/ (l - e^/(5o,fc-i) 



sup 



-log77Q"(l)- log 7700(G) 
n 



< b/n 



(4.18) 



for some finite constant b < 00. For a more thorough discussion on the stability 
properties of the semigroup ^>o,ra and the limiting measures ?7oo! we refer the reader 
to [Ij. Our next objective is to transfer these stability properties to the one of the 
sequence r/„. First, using ()4.18p . we readily prove that 

lim ilog7„(l) = logr?oo(G') 

n— >oo fi 

Next, we simplify the notation and we set a„ := Qo,n (7o(l);%) and 

Cn '■ — co rt- Using 

(j4.1ip . we find that for any n > 1 



a ^ \\r]n - r/oo||tv < e ^" + (1 - q„) ^ 
Recalling that 
we also obtain that 

1 



0<p<n So<<7<n 



E 

0<p<n 



El 



<q<n 



0<q<n 



1/r 



E 

0<p<n 



CriG 



-Xpr 



1/r 



< 



Eo<a<n 9- 



l/r 



0<p<n 



1/r 



(4.19) 



for any r > 1. We conclude that 



1 , 



0<p<n Eo<g<n 



and therefore 



< e-^- + /(I - e-'^g+)'/^ ^„^oo 



4.3 Stability and Lipschitz regularity properties 

We describe in this section a framework that allows to transfer the regularity properties of 
the Feynman-Kac semigroups <l>p,„ introduced in ()3.2p to the ones of the semigroup Tp^n 
of the sequence (7^(1), ^/n)- Before proceeding we recall a lemma that provides some weak 
Lipschitz type inequalities for the Feynman-Kac semigroup ^p,n in terms of the Dobrushin 
contraction coefficient associated with the Markov transitions Pp^n introduced in (|3.2p . The 
details of the proof of this result can be found in [4] or in [5] (see Lemma 4.4. in |5], or 
proposition 4.3.7 on page 146 in [3]). 
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Lemma 4.4 ([5]) For any < p < n, any r], ^ £ V{Ep) and any f £ Osci(-E„), we have 

li'^pAf^) - '^pAv)] if)\ < 2 qln l3iPp,n) |(/u - v)T^P,nM)\ (4-20) 

for a collection of functions 'Dp^n,r]{f) £ Osci(£'p) whose values only depend on the parame- 
ters (p, n, T]). 

Proposition 4.5 For any < p < n, any r],r]' G V{Ep) and any f G Osci(£'„), there exits 
a collection of functions T>p^n,r)'{f) € Osci(£^p) whose values only depend on the parameters 
{p, n, rf) and such that, for any m € Ip, we have 

|[r2^„(m,r?)-r2^„(m,V)] (/)| 

(4.21) 

< 2 a^^„ qp^n [qp,n l3{Pp,n) \{V - V')'^p,n,r,'if) \ + /3p,n | (?? " V)/lp,n,r?' |] 

with the collection of functions hp^n,'q> = }7%'"''(i) ^ C)sci(£'p) and the sequence of 

parameters Cp^n o-i^d /3p „ defined below 

<n := «p,n("^+(p)) and /3p,„ := ^ ^^-^ P{Pq,n) (4.22) 

p<g<n ^P<r<n Cr,n 

Before getting into the details of the proof of proposition 14.51 we illustrate some conse- 
quences of these weak functional inequalities for time- homogeneous models (£'„, G^, M„, ^„) = 
{E,G,AI,fj,) in the three scenarios discussed in (|4.13p . (|4.14p . and (I4.16p . 

1. When G{x) = 1 for any x £ E, we have 

^P,n(r/) = = 1/2 Cp,„ = /i(l) qp,n = l "p,™ < 1 

Let us assume that there exist a < oo and < A < oo such that /3(M") < ae""**"" for 
any n > 0. In this situation, we prove using (I4.2ip that 

|[r2 „(m,77) -r2 Jm,r?')] (/)| < 2ae-^("-P) |(^ - r?)Pp,„,,,(/)| 

2. When < 1 and when the mixing condition (M)^ stated in (j4.ip is satisfied for some 
A; and some fixed parameters ep = e, we have seen in (j4.15p that 

sup Q;„(m) < 1 A with d = ((7o(l)//i(l)) V (1 - g+y^) ^c^e"' 

Furthermore, using the estimates given in (]4.3p and (j4.4p . we also have that 

qp,n<Sk/e /3p,„ < 1 and /3(Pp,n) < a e~-^ with (a, A) given in gUD 

In this situation, we prove using ()4.2ip that 

\[rlArn,v)-rlnim,v')] if)\ 

< 2 [l A (d 5?"''^)] (<5fc/e) [(5fc/e) a e-^("-f) - r?)Pp,„,^'(/) | + |(^ " ^)^,nY|] 
Notice that for {n — p) > log (d) / log (1/(7+), this yields 
|[r2jm,7?)-r2 „(m,r?')] (/)| 

< ao e-^o("-P) Km - 7?)Pp,„,v(/)| + ai e-^i("-f) |(m - 7?)/ip,„,v| 

with 

ao = 2ad{5k/ef ai = 2d{5k/e) Aq = A + log (l/c/+) and Ai = log (1/5-+) 
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3. When g_ > 1 and when the mixing condition {M)k presented in (j4.ip is met for some 
k and some fixed parameters ep = e > 0, then we use the fact that 

ap,„ < 1 qp,n < h/e and /3(Pp,„) < a e-^^"-^) with (a, A) given in 

Arguing as in (|4.19p . we prove that for any > ^ log 5+ 

from which we conclude that 

|[r2jm,,7)-r2 jm,r?0] (/)| 



with 



ao = 2a{6k/ef = 2gL{Sk/e)/{l - e-^^g+)'/'' Aq = A and Ai = log(5_) 



Now, we come to the proof of proposition 14.51 
Proof of proposition 14. 5t 

First, we observe that 

^l,nim,r]) - r2_„(m',7y') 

ap,n (m, 7?) [$p,n(r?) - Ep<q<n '^J'^lc^,^ ^lA'P, 



•p<Cr<.n 



Using the fohowing decomposition 

ah - a'h' = a'{b - b') + (a - a')b' + (a - a'){b - b') 
which is vahd for any a, a' , 6, b' G M, we prove that 
rp,„(m,??) -Tlnim',!]') 



+ 



+ [ap,n ("1, r]) - Up^n (m', rj')] - ^p,n{^')] 

For 171 = 171', using (j4.24p we find that 

= Op^n {m, rj) [^p,niv) - ^p,n{^')] 



+ 



We also notice that 



ap,n {m,rj) = - 



+ /ip,„/[mryQp,„(l)] 
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(4.23) 



(4.24) 



from which we easily prove that 

otp,n im,r]) - Op^n {m', ry') 

and therefore 

ap,n {m, T]) - ap^n {m, t]') = (ap,„ (m, r/') (1 - Qp,„ (m, ??))) [t] - rj'] ( ^'"'^n^ 

The proof of ap^„ {rn,rj) < ctp^ni^) is elementary. Prom the above decomposition, we prove 
the following upper bounds 

\ap,n {m, rj) - Op^n {m, rj') \ < ap.„(m) [rj - r]'] ( 



and 



|[r2^Jm,r?)-r2_Jm,V)] (/)| 
<<nM[|[^P,n(r?)-$p,„(r?')] (/)| 



+ 

This yields 



|[r2 „(m,ry)-r2 „(m,r?')] (/)| 

< a^,„(m) [|[<I>p,„(ry) - ^p^v')] (/)l + /5p,n |[r? - v'] {^^) 
The last formula comes from the fact that 

PiPq,n):= sup ||^>g,nH -$q,n(l^')||tv 



The proof of this result can be found in [4J (proposition 4.3.1 on page 134). The end of 
the proof is now a direct consequence of lemma l4^ This ends the proof of the proposition. ■ 



5 Mean field particle approximations 

5.1 McKean particle interpretations 

In proposition 14. H the evolution equation (j4.5p of the sequence of probability measures 
iln Vn+i is a combination of an updating type transition ry„ '^G„iVn) and an integral 
transformation w.r.t. a Markov transition M„^x,(7„(i),»7„) that depends on the current total 
mass 7n(l) and the current probability distribution The operator M„_j_i_(^^(i) j^^) defined 
in ()4.6p is a mixture of the Markov transition M^+i and the spontaneous birth normalized 
measure ]In+i- We let Sn,r]„ be any Markov transition from En into itself satisfying 
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The choice of these transitions is not unique. We can choose for instance one of the col- 
lection of transitions presented in (j2.10p . (j2.1ip and (j2.12p . Further examples of McKean 
acceptance-rejection type transitions can also be found in section 2.5.3 in [4j. By construc- 
tion, we have the recursive formula 

r]n+i = r]nKn+i,{'r„{i),,j„) with i^„+i,(^„(i),^„) = 5'„,^„M„+i^(^^(i) .^^) (5.1) 

with the auxiliary total mass evolution equation 

7„+i(l) = 7„(1) r/„(G„) + /i„+i(l) (5.2) 

As already mentioned in section [21 the sequence of probability distributions rjn can be 
interpreted as the distributions of the states Xn of a process defined, conditional upon 
(7n(l),r7n), by the elementary transitions 

P {Xn+i G dx I Xn) = i^„,(^„(i),^„) {Xn, dx) with r/„, = Law(X„) 

Next, we define the mean field particle interpretations of the sequence (7^(1), ?7n) given in 
()5.ip and (j5.2p . First, mimicking formula ()5.2p we set 

7^+1(1):= 7^(1) r?^(Gn) + Mn+i(l) and 7;^(/) = 7^(1) x r,^(/) 
for any / G B{En), with the initial measure '^^ = 70. It is important to notice that 

n 

7^(l)=7o(l) n <(G,) + J;M1) n <(G',)^7^(l)e/n 

0<g<n p=l p<q<n 

The mean field particle interpretation of the nonlinear measure valued model ()5.ip is an - 
valued process ^„ with elementary transitions defined in (j2.9p and (j5.ip . By construction, 
the particle evolution is a simple combination of a selection and a mutation genetic type 
transition 

~^ in = iCn)l<i<N S,n+1 

During the selection transitions ^„ each particle -w evolves according to the 

selection type transition ,^jv (^^, dx). During the mutation stage, each of the selected 
particles ~^ Cn+i evolves according to the transition 

M„+i (^Ar(i) ,jjv)(x,(i?/) := an {in Mn+i{x,dy) + (l - a„ (7,^^(1), r?^)) 7^„+i(dy) 

5.2 Asymptotic behavior 

This section is mainly concerned with the proof of theorem 13.21 In section 15. 2. H we discuss 
the unibiasedness property of the particle measures 7^ and their convergence properties 
towards 7„, as the number of particles N tends to infinity. We mention that the proof of the 
non asymptotic variance estimates (j3.1ip is simpler than the one provided in a recent article 
by the second author with F. Cerou and A. Guyader [3]. Section 15.2.21 is concerned with 
the convergence and the fluctuations of the occupation measures r/^ around their limiting 
measures rjn- 
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5.2.1 Intensity measures 



We start this section with a simple unbiasedness property. Recall that J-p^'^ stands for the 
(T-field generated by the random sequence {s!j^'^)a<k<p- 



Proposition 5.1 For any < p < n, and any f £ B{En), we have 



p<q<n+l 



(5.3) 



In particular, we have the unbiasedness property: E (7^ (/)) = ^n{f)- 
Proof: 

By construction of the particle model, for any / G B{En) we have 

E [ri^+iU) I 4^)) = r?^i^n+i,(,-(i),,-)(/) = Tl_,, {^^{lU^) (/) 

with the second component r^_|_]^ of the transformation r„_|_i introduced in 14.71 Using the 
fact that 

p2 IN,-, . N\ ( .^ ^ 7^ (1) Vn (Qn+l(/)) + ^n+l(/) ^ 7^(Qn+l (/)) + (/) 
u+lKlnK ),Vn)U) ^^A.(l)^A.(g„^,(l))+^„^^(l) 7A^(Q„+i (1)) + (1) 

and 

7^+l(l) = 7^(1) (Gn) + /^n+l(l) = 7;:'(Qn+l(l)) + /^n+l(l) 

we prove that 

IE (7^+i(/) I ) = E (7„^+i(l) v^^iif) I 4"^)) = 7n";i(l) E (v^^Af) I 4"^)) 

= ln{Qn+l{f))+^^n+l{f) 

This also implies that 

e(7;:;i(/)| 4"^!) =iE(7„^(Qn+i(/))| 4"^!) 

= 7^-l(QnQn+l(/)) + fln{Qn+l{f)) + l^n+l{f) 

Iterating the argument one proves (j5.3p . The end of the proof is now clear. ■ 

The next theorem provides a key martingale decomposition and a rather crude non 
asymptotic variance estimate. 

Theorem 5.2 For any n > and any function f G B{En), we have the decomposition 



N [7/r -7n] (/) = KiQpAf)) 

p=0 



(5.4) 



In addition, if the mixing condition (M)fc presented in (fJTl) is met for some k > 1 and some 
constant parameters Cp = e > 0, then we have for any N > 1 and any n > 1 



E 



7^(1) 

7n(l) 



< 



n + l 6l 



51 



— 1 + 

e2 e'^{N -I) 



n~l 



(5.5) 
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Before presenting the proof of this theorem, we would hke to make a couple of comments. 
On the one hand, we observe that the unbiasedness property follows directly from the 
decomposition (j5.4p . On the other hand, using Kintchine's inequality, for any r > 1, p > 1, 
and any / G Osci(-En) we have the almost sure estimates 



(N) 



p-1 



A detailed proof of these estimates can be found in [3], see also lemma 7.2 in [Ij for 
a simpler proof by induction on the parameter A^. From this elementary observation, and 
recalling that 7^ (1) £ for any n > 0, we find that 



1 

iVE(|[7^-7n] {f)\y <ar h 



for some finite constant 6„ whose values only depend on the time parameter n. 
Now, we present the proof of theorem 15. 2[ 
Proof of theorem 15. 2t 
We use the decomposition: 



By ([O]), we find that 

7^+i(/)-e(7;:;i(/) =^'?Vi(/)- [7^(^3"+i(/)) + Mn+i(/)] 



Since we have 



7^ (Qn+l(l)) + /in+l(l) = 7n (Gn) + /in+l(l) 



N, 



this implies that 

7^+i(/) - [7^(Q„+i(/)) +/.„+!(/)] = 7;:;i(i) 



7^(1) r?^(G„) + A.n+i(l) = 7^:^1(1) 

K(Qn+l(/))+Mn+l(/)] 



[7^(g„+l(l))+Mn+l(l)] 
= 7,^^+1(1) [r]n+lif) - r]n ^n+l,(7,^'(l),r,^)(/) 

and therefore 

7;^+i(/)-e(7^+i(/)| 4^)) =7/y+i(l) [^^+i(/)-r?^i^„+i,(,-(i),,-)(/)" 
Finally, we observe that 

IE (7^+l(/) I 4^)) - 7n+l(/) = l^{Qn+l{f)) - 7n(gn+l(/)) 

from which we find the recursive formula 

[7/y+l - 7n+l] (/) = 7^+1(1) - r?^i^n+l,(,-(l),,-)] (/) + [in " 7n] (Qn+l(/)) 

The end of the proof of ()5.4p is now obtained by a simple induction on the parameter n. 
Now, we come to the proof of (|5.5|) . Using the fact that 



E (7p^(lX(/^)) 7,^(lX(/(^^)) = IE (7p^(l)7f (lX(/(^^) IE (<(/^)) I -Fjl,)) 

= 
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for any < p < q < n, and any /(^) G B{Ep), and /^^^ G B{Eq), we prove that 
iV E - 7n(l)]') = (7p^(l)' E 



p=0 

Notice that 

1 ^ 1 1 / 7p(Qp,n(l)) 

7n(l)2 ^^(1)2 r?,(gp,„(l))2 1^ 7„(1) 

The r.h.s. estimate comes from the fact that 

7p(Qp,n(l)) _ 7p(l) '?p('3p,n(l)) 



7n(l) 7p(l) ??p(Qp,n(l)) + Ep<g<n lJ^qQq,n{^) 

Using the above decompositions, we readily prove that 



ap,n (7p(l),^p) < ap,n(7p(l)) 



N E 
with 
We set 



7^(1) ^ 

7n(l) 



<E<n(7p(l))'E 
p=0 



7p^(l) 

7.(1) 



E«(Qp,„(i))Vp-i) 



Qp,n(l) = <3p,n(l)/^p(Qp,n(l)) < 9p,i 



C/i^ := E 



7^ 

7n(l) 



then we find that N < a„ + ^ ^p,n 

p=0 



AT 
P 



with the parameters 



•= E (^P,n«P,n(7p(l))^ and 6p,„ := (gp,„a* „(7p(l))^ 

p=0 

Using the fact that bn,n < 1, we prove the following recursive equation 

1 1 AT 

bp 

.n 



<<+Y. <n < With 

0<p<n 



N-1 ^^'--■=N-1 



Using an elementary proof by induction on the time horizon n, we prove the following 
inequality: 





n 








E< E ^"(^) 


+ 


E ^''('^) 




_P=1 ee{p,n} 




_ee{0,n) 



In the above display, {p,n) stands for the set of all integer valued paths e = (e(/))o</<jt of 
a given length k from p to n 



eo 



= p < ei < . . . < ejfc_i < efe = n and fe^(e) = JJ 6, 



AT 

e(i-l),e(0 



KKfe 



We have also used the convention 6^(0) = Yl^ = 1 and (n, n) = {0}, for p = n. Recalling 
that 7^ = 70; we conclude that 

n 

f^n^<E< E b'^ie) 
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We further assume that the mixing condition (M)jt presented in (|4.ip is met for some 
parameters A; > 1, and some constant parameters = e > 0. In this case, we use the fact 
that 

0/-p^n — ^ and Qp^n 

to prove that 

sup <{n + 1) {6k/ef/{N - 1) and sup < {6k/ef/{N - 1) 

0<p<n 0<p<n 

Using these rather crude estimates, we find that 

0<p<n 1=1 ^ / \ \ I J 

and therefore 

(N-l) £2 V e^N -1)) 
This ends the proof of the theorem. ■ 




5.2.2 Probability distributions 

This section is mainly concerned with the proof of the Lj-mean error estimates stated in 
(|3.9p . We use the decomposition 

(7;r(l),^^) - (7n(l),r/„) = [ro,„ (70^(1),%^) -ro,„ (70(1),%)] 

n 

+ E (7p^(l).0 - rp-i,n (7p^-i(l),^p^-i)] (5.7) 
p=i 

to prove that 

A'' 

= [r2„ (70^(1),%^) - rl^ (70(1),%)] + e;=i [r^,„ (7p^(i),0 - r^_i,„ (7,^.1(1), r?f_0] 

Using the fact that 

rp_i,„(m, 77) = Tp^n {Tp{m, rj)) r^_i_„(m, r/) = T^ „ (rp(m, rj)) 
we readily check that 

(7p^-i(l),^p^-i) = (7p''-i(l)r/p^-i(G'p-i)+;Up(l),vI/G,., (!)<_,)) 

= (7^(1)' ^^-l^p,(7,^_,(l),<„i)) 
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Since we have 7(^(1) = Aio(l) =7o(l); one concludes that 

v!!-vn= [rg,, (70(1),%^) - (70(1),%)] 

n 
p=l 

Using the fact that 7^ (1) E Ip, for any p > 0, the end of the proof is a direct consequence of 
lemma 14.51 and Kintchine inequality. The proof of the uniform convergence estimates stated 
in the end of theorem 13.21 are a more or less direct consequence of the functional inequalities 
derived at the end of section [5]3j The end of the proof of the theorem 13.21 is now completed. 

We end this section with the fluctuations properties of the A^-particle approximation 
measures 7^ and rj^ around their limiting values. Using the type of arguments as those 
used in the proof of the functional central limit theorem, theorem 3.3 in [6], we can prove 
that the sequence {W^)n>o defined in ()3.12p converges in law, as tends to infinity, to the 
sequence of n independent, Gaussian and centered random fields (M^n)ri,>o with a covariance 
function given in (j3.13p . Using the decompositions (15. 4p and 

- ,„(/) = ^ ([7." - T,J (^(/ - ,.(/») ) 
by the continuous mapping theorem, we deduce the functional central limit theorem 13.31 



6 Particle approximations of spontaneous birth measures 

Assume that the spontaneous birth measures are chosen so that /i„ <C A„ for some 
reference probability measures A^^ and that the Radon Nikodim derivatives Hn — d^jij dX^i 
are bounded. For any n > 0, we let A^' := '}2f=i ^Cn empirical measure associated 

with N' independent and identically distributed random variables (Cn) i<j<Ar with common 

distribution A^. We also denote by fi^' the particle spontaneous birth measures defined 
below 

Vn > {dXn) ■■= Hn{Xn) Xn i^Xn) 

In this notation, the initial distribution ryo ^ind the initial mass 70 are approximated by the 
weighted occupation measure r/^ := ^Ho{^q ) and 7,^ (1) := A^ {Hq). 

We let 7^ and rj^ the random measures defined as 7n and rjn by replacing in (12. 2[) the 
measures Hn by the random measures fi^ , for any n > 0; that is, we have that 

1^' =ltlQn + l^''n' and e'(/n) = 7n^'(/n)/7r(l) 

for any G B{En)- By construction, using the same arguments as the ones we used in the 
proof of (j4.8p , we have 

~N' _ \ " ,,Af'n 

7n — 2-^ ^P'" 
0<p<n 

This yields for any / £ 13 (En) the decomposition 



~N' 

In -In 



(/)= E Yp' -i^)^QpM)= E [Af - Ap] (i^p Qp,„(/)) 



0<p<n 0<p<n 



Several estimates can be derived from these formulae, including Lp-mean error bounds, 
functional central limit theorems, empirical process convergence, as well as sharp exponential 
concentration inequalities. For instance, we have the unbiasedness property 



E(7f(/))=7n(/) 
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and the variance estimate 



N E 



7r(/)-7n(/) 



K {HpQp,n{f) - \{HpQp,n{f))Y 

0<p<n 



Using the same arguments as the ones we used in (j5.6p , we prove the following rather crude 
upper bound 



N E 



7n(l) 



Vnif) 



0<p<n, 



We illustrate these variance estimates for time homogeneous models {En, Gn,Hn, Mn, fJ-n) = 
{E,G,H,M,^), in the three situations discussed in (j4.13p . (I4.14|) . and (I4.16j) . We further 
assume that the mixing condition {M)^ presented in (|4.ip is met for some parameters A; > 1, 
and some e > 0. In this case, we use the fact that qp^n < (^fc/e, to prove that 



N E 



7n(l) 



rinU) 



E [<n(7p(l))/7p(l)]' 
0<p<n 



with some constant c := {\\H\\ /i(l) ((5fc/e)^). 

1. When G{x) = 1 for any x £ E, we have 7p(l) = 7o(l) + /^(l) P- Recalling that 
Q^p,n(7p(l)) < 1, we prove the uniform estimates 



supE 

n>0 



7n(l) 



Vnif) 



<cJ](7o(l) + Ml)p)" 

p>0 



2. When 5+ < 1 and when the mixing condition {M)^ stated in (j4.ip is satisfied, we have 
seen in ()4.15p that 

<n(7p(l)) < 1 A (di 5?"^^) and inf 7„(1) > da 

for some finite constants di < 00 and ^2 > 0. From previous calculations, we prove 
the following uniform variance estimates 



N supE 

n>0 



Jim 

7n(l) 



^n(/) 



< 



2 „2p 



p>0 



3. When (/_ > 1 we have seen in (j4.16p that 7n(l) > d g^l for any n > no, for some finite 
constant d < 00 and some no > 1 so 



N supE 

n>0 



TfCf) 

7n(l) 

















J 





0<p<no 



n>no 
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